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Abstract. We use a space-based plume height climatology derived from observations 
made by the Multi-angle Imaging SpectroRadiometer (MISR) instrument aboard the NASA 
Terra satellite to evaluate the ability of a plume-rise model currently embedded in sev- 
eral atmospheric chemical transport models (CTMs) to produce accurate smoke injec- 
tion heights. We initialize the plume-rise model with assimilated meteorological fields from 
the NASA Goddard Earth Observing System and estimated fuel moisture content at the 
location and time of the MISR measurements. Fire properties that drive the plume-rise 
model are difficult to estimate and we test the model with four estimates for active fire 
area and four for total heat flux, obtained using empirical data and Moderate Resolu- 
tion Imaging Spectroradiometer (MODIS) fire radiative power (FRP) thermal anoma- 
lies available for each MISR plume. We show that the model is not able to reproduce 
the plume heights observed by MISR over the range of conditions studied (maximum 
r 2 obtained in all configurations is 0.3). The model also fails to determine which plumes 
are in the free troposphere (according to MISR), key information needed for atmospheric 
models to simulate properly smoke dispersion. We conclude that embedding a plume- 
rise model using currently available fire constraints in large-scale atmospheric studies re- 
mains a difficult proposition. However, we demonstrate the degree to which the fire dy- 
namical heat flux (related to active fire area and sensible heat flux), and atmospheric 
stability structure influence plume rise, although other factors less well constrained ( e.g ., 
entrainment) may also be significant. Using atmospheric stability conditions, MODIS FRP, 
and MISR plume heights, we offer some constraints on the main physical factors that 
drive smoke plume rise. We find that smoke plumes reaching high altitudes are charac- 
terized by higher FRP and weaker atmospheric stability conditions than those at low 
altitude, which tend to remain confined below the BL, consistent with earlier results. We 
propose two simplified parameterizations for computing injection heights for fires in CTMs 
and discuss current challenges to representing plume injection heights in large scale at- 
mospheric models. 


1. Introduction 

Biomass burning is a significant source of trace gases and 
aerosols to the atmosphere, with the potential to alter its 
chemical and radiation properties over extensive regions. 
Fires are also important sources of heat, which generates 
strong updrafts above the fire. The energy emitted by fires, 
in combination with atmospheric conditions, determines the 
vertical distribution of the fire emissions in the atmosphere 
near the fire source (i.e., injection height). Global and re- 
gional chemical transport models (CTMs) have been used to 
study the impact of fire emissions on air quality and atmo- 
spheric composition [e.g. Pfister et al., 2006 ; Turquety et al., 
2007; Spracklen et al., 2007]. Thus, knowledge of the injec- 
tion heights is fundamental to accurately representing fire 
emissions in CTMs. However, representing injection heights 
in models is a difficult task, complicated by the sparseness 
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of plume height measurements to validate the calculated 
values. The space-based plume height climatology derived 
from observations made by the Multi-angle Imaging Spec- 
troRadiometer (MISR) instrument aboard the NASA Terra 
satellite provides a unique dataset for extensive validation 
of injection heights. We use these data to evaluate a one di- 
mensional (1-D) plume-rise parameterization currently used 
to simulate injection heights in CTMs, and to investigate 
the main physical factors that determine smoke plume rise. 

Fire emissions reach a wide range of altitudes in the at- 
mosphere. Plume heights have commonly been measured for 
isolated fires on a case-by-case basis. For example, smoke 
plumes have been observed completely confined within the 
atmospheric boundary layer (BL) [e.g. Trentmann et al, 
2002] and at different altitudes in the free troposphere (FT) 
(~2 to 7 km) [e.g. Wofsy et al, 1992; de Gouw et al, 
2006]. Wildfire plumes have also been detected within 
the lower stratosphere (>10 km) [e.g. Fromm et al., 2005; 
Damoah et al., 2006; Dirksen et al., 2009], although this 
phenomenon occurs very sporadically and is associated with 
pyro-convective events triggered by intense heat from the fire 
and very unstable atmospheric conditions. 

Recent advances in space-based observations have pro- 
vided a tool for identifying smoke injection heights at 
regional-to-global scales. Smoke plume heights have been 
retrieved using space-borne lidar observations [e.g. Labonne 
et al., 2007; Raffuse et al., 2012], stereo imaging [e.g. Kahn 
et al, 2007, 2008; Val Martin et al., 2010; Mims et al., 
2010 ; Tosca et al., 2011] and aerosol index (Al) measure- 
ments [Guan et al., 2010]. Labonne et al. [2007] studied 
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hundreds of plumes distributed around the world using data 
from the Cloud-Aerosol Lidar with Orthogonal Polariza- 
tion (CALIOP) on board the Cloud- Aerosol Lidar Infrared 
Pathfinder Satellite Observation (CALIPSO) satellite and 
found that smoke is mostly confined within the BL. In con- 
trast, Kahn et al. [2008] found that, based on stereo-height 
retrievals from MISR, about 5-18% of the plumes reached 
the FT over Alaska and the Yukon territories in 2004 These 
results were supported by the comprehensive study of Val 
Martin et al. [2010], who analyzed about 3000 plumes across 
North America during 2002 and 2004-2007, and found that 
an important fraction (4-12%) of plumes inject smoke into 
the FT over North America. Following the case study anal- 
ysis of Kahn et al. [2007], Val Martin et al. [2010] showed 
that smoke from the plumes reaching the FT typically get 
trapped within atmospheric stable layers, when they exist; 
otherwise smoke tends to spread-out vertically. A MISR 8- 
year plume height climatology for tropical forest and peat- 
land fires over Borneo and Sumatra found that less than 4% 
of plumes reached the FT, suggesting that the direct injec- 
tion of smoke into the free troposphere is not an important 
mechanism for vertical mixing of emissions over equatorial 
Asia [Tosca et al, 2011]. 

In addition to MISR smoke plumes, Val Martin et al. 
[2010] also analyzed smoke clouds retrieved from MISR. 
Smoke clouds are smoke without a detectable fire origin and 
constitute a later stage of smoke plume evolution. This anal- 
ysis showed that about 35% of smoke clouds reach the FT, a 
larger fraction than those of smoke plumes. The authors sug- 
gested that smoke from plumes may reach higher altitudes 
later in the day in part because fire intensity increases in the 
late afternoon, and in part because of atmospheric processes 
unrelated to the fire itself, such as convection and/or advec- 
tion. High altitude plumes have also been detected using an 
empirical method based on Al from the Ozone Monitoring 
Instrument (OMI) and the Total Ozone Mapping Spectrom- 
eter (TOMS). Guan et al. [2010] identified 180 plumes above 
8 km around the globe for a 10-year period. These findings 
are consistent with the previous documented detections of 
smoke in the free troposphere and lower stratosphere from 
isolated cases [e.g. Fromm et al., 2005; Damoah et al, 2006; 
Dirksen et al, 2009]. However, it is important to note that 
these two studies are biased to dense high altitude smoke 
plumes, as dense smoke is easier to visualize by the MISR 
operator [ Val Martin et al., 2010] and to detect with the Al 
analysis [Guan et al, 2010]. 

Smoke injection height is an important input in CTMs, 
as it has a strong influence on smoke dispersion [e.g. Colarco 
et al, 2004]- Smoke lofted into the free troposphere is often 
transported hundreds or thousands of kilometers downwind 
because of the higher wind speeds at higher altitudes [e.g. 
Val Martin et al., 2006; Dirksen et al., 2009], whereas smoke 
trapped within the boundary layer typically is well mixed 
and remains near the source region [e.g. Trentmann et al., 
2002]. Removal processes are also more efficient in the BL, 
whereas pollutants’ lifetime increases when they are trans- 
ported in the free troposphere [Chatfield and Delany, 1990]. 
At present, a wide range of arbitrary procedures is used to 
represent the vertical distribution of fire emissions in CTMs. 
Typically, a constant emission height is used. Emissions are 
either released initially within the BL [e.g. Lamarque et al, 
2003; Colarco et al., 2004], released at specified altitudes 
based on a simple empirical height-fire intensity relation- 
ship [e.g. Lavoue et al, 2000; Wang et al., 2006], uniformly 
mixed throughout the troposphere [e.g. Pfister et al., 2006], 
or a pre-selected fraction is released within and above the 
BL [e.g. Leung et al., 2007; Turquety et al., 2007; Generoso 
et al, 2007; Hyer et al., 2007]. However, all these schemes 
neglect the effects of both the buoyancy generated by the 
fire and the atmospheric conditions, with the potential of 
introducing significant errors in CTM simulations. 

Recently, more elaborate plume-rise parameterizations 
based on empirical and physical approaches have been in- 
corporated into CTMs to take into account fire buoyancy 


and atmospheric conditions. Table 1 lists some of these 
parameterizations and the corresponding embedded CTM. 
For example, DAYSMOKE and Briggs are two empirical- 
statistical plume rise schemes [Achtemeier et al, 2011; 
Briggs, 1969, respectively] embedded in the regional Com- 
munity Multiscale Air Quality (CMAQ) model [e.g. Liu 
et al., 2010; Achtemeier et al, 2011; Raffuse et al, 2012]. 
An evaluation of DAYSMOKE with six prescribed fires 
showed that the simulated injection heights were near the 
lower end of the range of observed plume heights. Simi- 
larly, Raffuse et al. [2012] found that simulated heights of 
smoke plume with the Briggs scheme were systematically 
lower than the values detected with MISR and CALIOP. 
The MISR-Briggs intercomparison also showed that mod- 
eled injection heights had much greater variability than the 
MISR heights, and the correlation between the two data sets 
was weak. These results suggest that accurately estimat- 
ing plume heights with empirical schemes is a big challenge, 
since they relay on static meteorological conditions (e.g., 
wind speed and atmospheric stability), uncertain calcula- 
tions of sensible heat and neglect the additional buoyancy 
generated by latent heat [e.g. Raffuse et al, 2012; Liu et al, 
2010; Stein et al., 2009]. 

One of the most widely used parameterizations is the 
plume-rise model developed by Freitas et al. [2007], and 
recently modified to include the effect of horizontal winds 
[Freitas et al., 2010]. The plume-rise model uses a prognos- 
tic scheme for simulating plume heights for tropical fires. 
The model determines fire buoyancy using the total sensi- 
ble heat flux and area of each fire, and takes into account 
the additional buoyancy from latent heat. The plume-rise 
model was first embedded in the regional CATT-BRAMS 
model [Freitas et al., 2007] and the global CAM model [Guan 
et al., 2008] for exploratory case studies in the tropics. Fre- 
itas et al. [2007] demonstrated, using Atmospheric Infrared 
Sounder (AIRS) and aircraft observations, that explicitly 
accounting for injection processes significantly improves re- 
gional atmospheric model predictions of tropospheric CO 
across the central Amazon. At present, the validation of 
this plume-rise model has focused only on three case stud- 
ies: the Quinault Fire described by Trentmann et al. [2002] 
[Freitas et al., 2007] and two simulated fires with the 3-D Ac- 
tive Tracer High Resolution Atmospheric Model (ATHAM) 
[Freitas et al., 2010]. 

More recently, the Freitas et al. [2007] plume-rise model 
has been embedded into WRF-Chem and used to study fires 
in North America [Pfister et al., 2011; Sessions et al, 2011; 
Grell et al., 2011], using the heat fluxes and typical fire sizes 
for vegetation types in the tropics from Freitas et al. [2006] . 
Sessions et al. [2011] compared injection heights simulated 
by the WRF-Chem online plume-rise model to MISR heights 
obtained during the ARCTAS campaign over Siberia and 
Canada in 2008 and showed a poor correlation (r 2 <0.2) 
between the model and MISR plume heights. The authors 
also evaluated the long-range transport simulated by WRF- 
Chem with three injection height schemes ( i.e ., plume-rise 
model, injected at 3-5 km and released within the BL). They 
found that the plume-rise model scheme produced the best 
match to AIRS observations, although the difference from 
results with the 3-5 km injection scheme was not significant. 
This underlines the complexities of representing plume in- 
jection heights in atmospheric models without compensating 
for errors within CTMs themselves, such as emissions and 
convective transport processes [e.g. Chen et al., 2009; Hyer 
and Reid, 2009]. 

In this work, we evaluate the plume-rise model described 
by Freitas et al. [2007, 2010] using the MISR plume height 
climatology presented in Val Martin et al. [2010]. We first 
develop a set of approaches to determine the fire character- 
istics needed to initialize the model (i.e., active fire area and 
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total heat flux). We also use the MISR plume heights, in 
combination with assimilated meteorological fields from the 
NASA Goddard Earth Observing System (GEOS) and Mod- 
erate Resolution Imaging Spectroradiometer (MODIS) fire 
radiative power thermal anomalies available for each MISR 
plume, to constrain smoke plume rise, independently of the 
plume-rise model. On the basis of these results, we recom- 
mend how to treat injection heights in CTMs and discuss 
current limitations and challenges to representing plume in- 
jection heights in large-scale atmospheric studies. 

2. Data and methods 

2.1. Aerosol smoke plume observations 

We use observations of near-source plume heights over 
North America for the 2002 and 2004-2007 fire seasons, ob- 
tained with the MISR instrument. The plume height dataset 
was developed with the MISR Interactive eXplorer (MINX) 
software package (http://www.openchannelfoundation.org/ 
projects/MINX). To retrieve smoke plume heights, MINX 
first locates fires in the MISR field of view using MODIS 
fire pixels and an operator then identifies the source and 
digitizes the boundaries of the plumes and determines 
the plume direction of transport. In addition to smoke 
plume heights, MINX also retrieves wind speeds and albe- 
dos, and provides aerosol properties (e.g., Angstrom expo- 
nent) from the MISR standard aerosol product, and to- 
tal fire radiative power (FRP) and brightness temperatures 
from the MODIS thermal anomaly product [Nelson et al., 
2008a]. Each smoke plume height retrieval is assigned a 
flag that describes the quality of the digitization output 
(i.e., ’’good”, ’’fair” or ’’poor”) [Nelson et al, 2008b]. The 
MINX plume data are publicly available at http://www- 
misr.jpl.nasa.gov/getData/accessData/MisrMinxPlumes/ in- 
dex.cfm. A detailed analysis of this climatology was pre- 
sented in Val Martin et al. [2010]. 

Plume-rise models typically provide a maximum injection 
height above the fire, as described in section 2.2. A number 
of statistics are used to describe the heights of MISR plumes, 
such as the median, mode and maximum of the vertical dis- 
tribution of the stereo-height retrievals [e.g. Kahn et al., 
2008; Val Martin et al, 2010; Tosca et al, 2011]. In this 
work, we use plumes with a data quality flag of ’’good” and 
modify the definition of plume height derived from MISR 
to be compatible with that simulated by 1-D plume-rise 
models. We define the MISR plume height as the maxi- 
mum of the stereo-height retrievals immediately above the 
fire source, whereas the heights in our previous work were 
aggregated over the entire plume. We take the fire source 
to be the area formed by the MODIS fire pixels clustered at 
the origin of the fire plume, and include a buffer around this 
cluster of 10% of the cluster radius. Following the quality 
protocol established for plumes digitized by MINX [Nelson 
et al., 2008a], we limit cloud contamination by setting the 
maximum height as the 99th percentile of the smoke pixels 
retrieved above the fire source. In addition, only plumes 
having more than five smoke retrievals above the fire source 
are considered. To characterize the fires associated with 
each plume, and to select cases with well-defined fires, we 
consider only those plumes having MODIS fire pixels clus- 
tered at the plume origin. We also exclude from the analysis 
pixels that are located more than 2 km away from the clus- 
ter. These additional constraints removed 62% of the good 
quality plumes in the dataset studied by Val Martin et al. 
[2010], leaving a total sample size of 835 plumes. 

Figures la and lb compare the vertical distributions of 
stereo-height pixels provided by the MISR dataset, and the 
stereo-height pixels that we consider to be above the fire 
source. The modification of the plume height definition 


does not significantly change the plume height analysis of 
Val Martin et al. [2010], it simply provides a more appropri- 
ate comparison between MISR-derived plume heights and 
model injection heights. For example, with the new defini- 
tion, the plume in Figure la has a best estimated maximum 
height above ground level of 3150 m and a total fire radiative 
power of 2100 MW, whereas these parameters are 3325 m 
and 2325 MW for the entire plume, used in the original 
dataset. For analysis purposes, a plume altitude of 3150 m 
above the fire source is more representative of a modeled in- 
jection height, as it does not take into account atmospheric 
processes unrelated to the fire itself that may affect the final 
plume rise; a total FRP of 2100 MW gives an estimate of the 
fire radiative heat output for a well-defined fire, and avoids 
including in the analysis fire pixels from areas burning near 
the plume origin, which the model cannot simulate. 

2.2. Plume-rise model 

We evaluate the 1-D plume-rise model described by Fre- 
itas et al. [2007] and Freitas et al. [2010]. This plume-rise pa- 
rameterization scheme was originally developed by Latham 
[1994], and assumes that smoke and trace gases are trans- 
ported upward due to positive buoyancy produced by the 
fire. The altitude reached by these emissions is a complex 
function of fire energy and local meteorology at the time 
of burning. Within the column above the fire there is a 
decrease in temperature due to radiative cooling, and effi- 
cient heat transport by convection. In addition, the plume 
dilutes and loses buoyancy due to entrainment of colder en- 
vironmental air into the smoke plume, and lateral entrain- 
ment due to the drag of the wind. An additional positive 
buoyancy is sometimes generated from latent heat, when 
water vapor condenses [Freitas et al, 2007]. The 1-D plume- 
rise model accounts for these processes, and is governed by 
equations based on energy conservation, vertical momentum 
conservation, and continuity of water phases [Freitas et al., 
2007, 2010]. The final height of the plume (i.e., the top of the 
plume or “injection height”) is determined as the height at 
which the vertical velocity of the air parcel is less than 1 m/s 
[Freitas et al., 2006]. Figures lc and Id show three output 
parameters estimated by the model (plume vertical velocity, 
vertical mass distribution and acceleration buoyancy) at the 
locations of two MISR plumes, one high-altitude plume ob- 
served in central Canada in 2005 and one low-altitude plume 
observed in Texas in 2006. In these cases, the model predicts 
that the vertical velocities of the simulated plumes diminish 
with altitude, and the plumes reach maximum heights above 
ground level (agl) of 3500 m and 850 m, respectively, similar 
to those observed by MISR for the actual plumes (3150 m 
and 596 m). 

2.3. Input data for plume-rise model 

The plume-rise model is driven by four main inputs: ac- 
tive fire area, total fire heat flux, fuel moisture content of the 
burning vegetation, and meteorological variables (i.e., rela- 
tive humidity (RH) , temperature and wind speed profiles), 
as discussed by Freitas et al. [2007]. The model calculates 
the buoyancy flux generated by the fire using the sensible 
heat flux, assumed to be 55% of the total fire heat flux, 
and the radius of the plume, which is estimated from the 
active fire size. Typically, the 1-D plume model is embed- 
ded in each column of a 3-D host model, in which the 3-D 
model feeds the plume model the environmental parameters 
[e.g. Freitas et al., 2007; Pfister et al., 2011; Sessions et al, 
2011]. In this work, we run the plume-rise model at the lo- 
cation of the MISR plumes, independently of a 3-D model. 
We initialize the model with assimilated GEOS meteoro- 
logical observations, and estimate the fuel moisture content 
using the Canadian Fire Weather Index (FWI) system [Van 
Wagner, 1987] at the location and time of the MISR plume 
observations. Additionally, we used several approaches to 
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determine the active fire area and the total heat flux for 
the MISR plume fires. These parameters are difficult to es- 
timate, and large uncertainties exist in the data currently 
available [e.g. Hyer and Reid, 2009; Schroeder et ai, 2010]. 
We discuss them separately in section 2.4. 

For meteorological conditions, following the work of Val 
Martin et al. [2010], we used the GEOS meteorological 
fields that were archived for running the GEOS-Chem global 
model, with degraded resolution of 2° latitude by 2.5° lon- 
gitude. The temporal resolution of the fields is 6 hours 
(3 hours for surface variables and mixing depths). We ex- 
tracted RH, temperature and wind speed data from GEOS-4 
for 2002 and 2004-2006, and GEOS-5 for 2007. GEOS-4 has 
55 vertical levels between the surface and 0.01 hPa (includ- 
ing 5 levels up to 2 km), whereas GEOS-5 has 72 vertical 
levels, with 14 levels up to 2 km. Due to the difference be- 
tween the horizontal resolution of GEOS and MISR (200 km 
versus 1.1 km), there was in some cases a mismatch between 
the terrain elevations. To avoid a bias in the results due 
to the terrain difference, we excluded the observations for 
which the difference between the MISR and GEOS terrain 
elevation exceeded 250 m. This resulted in the removal of 
about 30% of the plume cases, leaving a total of 584 plumes. 
In addition, to assess the relationship between the injection 
height given by the model and the atmospheric structure, 
we used the atmospheric boundary layer heights provided 
with the GEOS fields [ Lucchesi , 2007]. 

Vegetation fuel moisture content at the fire locations was 
estimated as the average of the fine fuel and duff mois- 
ture content, obtained from the Fine Fuel Moisture Code 
(FFMC) and the Duff Moisture Code (DMC), respectively, 
and calculated using surface meteorological observations 
with the FWI system. Typical fuel moisture content at the 
location of the MISR plumes was 40.3% (reported median 
of all observations), which is similar to the 40% employed 
by Trentmann et al. [2006] for an extra-tropical forest fire 
simulation. The fuel moisture values used here are some- 
what larger than the 10% used as default in the plume-rise 
model. However, we found that the injection height results 
are not particularly sensitive to the fuel moisture content. 

2.4. Fire buoyancy configurations 

To evaluate the performance of the plume rise model rel- 
ative to the fire buoyancy sources, we tested four approaches 
for determining the total heat flux and four for the active 
fire area. Table 2 summarizes the 16 possible configurations. 

2.4.1. Total Heat flux 

The approaches for total heat flux from the fires are 
named FT06, FUEL, Dual-Band and FRPxlO, as shown in 
Table 2. A summary of the heat flux values estimated using 
these approaches is presented in Table 3. 

Freitas et al. [2006] adopted upper and lower limit heat 
flux values for three types of vegetation in Brazil (grassland, 
woody savanna, and tropical forest), based on the fire liter- 
ature for that region. For our FT06 approach, we classify 
vegetation in three groups: grassland, shrubland and for- 
est, and use the values in Freitas et al. [2006], as these were 
adopted in previous modeling work for the extratropics [e.g. 
Pfister et al., 2011; Sessions et al, 2011; Grell et al., 2011]. 
We run the plume-rise model with the upper and lower limit 
values and report the average of the results, with the excep- 
tion of grassland, for which only one value is given in Freitas 
et al. [2006]. 

The FRP approach is based on the fire radiative power 
reported for each fire pixel within the MISR plume dataset. 
Previous studies determined that fire radiative energy is 
about 10-20% of the total fire heat energy, based on the 
heat content of the fuels, and measurements of fire bright- 
ness temperature and spectral ernittance [Wooster et al., 


2005; Freeborn et al, 2008]. We thus consider the total 
heat flux released in each fire as the total FRP measured 
in each plume multiplied by 10. It is important to note that 
the MISR plume database provides the FRP values in MW. 
However, we consider the MODIS FRP as MW per 1 km 
pixel (thus in units of W/m 2 ). This assumption is appro- 
priate in our analysis because the narrower MISR swath lies 
in the center of the MODIS/Terra swath. Thus, MODIS 
pixels are distorted very little near the center where MISR 
observes. The maximum correction to adjust the area to the 
scan angle applied to the MODIS pixels is 15% and this an- 
gle correction is used in a small fraction of the data [David 
Nelson, JPL, personal communication]. 

The Dual-Band approach is based on an implementa- 
tion of the Dozier [1981] bi-spectral algorithm applied to 
the original MODIS Level IB calibrated radiance data 
(MOD021km). Fire pixels within the MOD021km gran- 
ules are first detected with an active fire detection algo- 
rithm that has been designed for use with the forthcoming 
Sentinel-3 Sea and Land Surface Temperature Radiometer 
(SLSTR), and whose performance has been recently evalu- 
ated using MODIS and ASTER data [Wooster et al, 2012]. 
Each detected fire pixel from the MOD021km granule is as- 
signed to a particular ’’fire”, based on the analysis of clus- 
ters of spatially contiguous fire pixels. This is an approach 
similar to that used previously with data from the BIRD 
Hotspot Recognition Sensor [Wooster et al., 2003; Zhukov 
et al, 2005], and designed to minimize some of the prob- 
lems of the dual band method, for example those related to 
inter-channel spatial misregistration effects [Zhukov et al., 

2005] . Each fire cluster is then subject to analysis using 
the Dozier [1981] ’’dualband” approach, using the specific 
method taken by Zhukov et al. [2005]. The output for each 
cluster is an estimate of the effective active fire temperature 
and active fire area, based on the fire clusters signals in the 
MODIS MIR and LWIR spectral bands, along with an esti- 
mate of the clusters total FRP. As in the FRP approach, we 
consider the total heat flux released in each fire as the total 
FRP multiplied by 10. 

Our last approach uses fuel consumption to estimate the 
heat flux. We use an improved fuel consumption map for 
western North America [Spracklen et al., 2009], and our own 
estimates for boreal North America. A description of the 
fuel consumption map is found in Appendix A. We refer to 
this approach as FUEL. We obtain heat flux values (i.e., 
W/m 2 ) from fuel consumption (i.e., kg/m 2 ) by assuming 
that the average burning time of the fire is two hours, and 
using a standard heat of combustion for extra-tropical veg- 
etation of 18.7e6 J/Kg [Trentmann et al., 2006]. 

2.4.2. Active fire area 

We apply four approaches to estimating active fire area 
for each plume, based on MODIS data available for each 
MISR plume. It is well recognized that a simple count of 
fire pixels at the nominal 1-kin MODIS resolution tends to 
vastly overestimate the active fire area [e.g. Giglio et al, 

2006] . Our goal is therefore to retrieve fire areas at the 
sub-pixel level. Table 3 summarizes the statistical parame- 
ters for active fire areas by vegetation type obtained in each 
approach. 

The first two approaches, WRF-Pxl and FLAMBE-Pxl, 
are based simply on the number of MODIS fire pixels, and 
differ only in the area assigned to each pixel. The WRF- 
Pxl approach clusters the pixels within each plume, assum- 
ing that the mean area burned within each fire pixel is 
0.5 km 2 . This assumption is used in the WRF-Chem online 
plume-rise configuration [Gabriele Pfister, NCAR, personal 
communication]. In WRF-Chem, fire pixels are aggregated 
within each model grid box [e.g., 12x12 km in Pfister et al. 
[2011]). The FLAMBE-Pxl configuration is based on the 
Fire Locating and Modeling of Biomass Burning Emissions 
(FLAMBE) dataset [Reid et al., 2009]. FLAMBE provides 
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carbon and aerosol emissions at hourly intervals, based on 
the WF-ABBA and MODIS active fire products. It sim- 
ulates diurnal variability of emissions and reports burned 
area taking into account the diurnal variability of fires. For 
a 1 km MODIS fire pixel, FLAMBE assumes 0.625 km is 
burned, and partitions it into 24 hourly segments (i.e., 0 km 
for nighttime hours, 0.05 km for daytime hours and 0.025 km 
for transitions hours) [Edward Hyer, NRL, personal com- 
munication]. We therefore consider that, at the time of the 
MISR observations (11:00-14:00 local time over North Amer- 
ica), 0.05 km on average is burning within each fire pixel. We 
then sum all fire pixels within each plume. The WRF-Pxl 
approach returns active fire areas that are exactly a factor 
of 10 larger than the FLAMBE-Pxl approach, as shown in 
Tabic 3. 

The third approach, scaled-FRP, is based on the FRP 
data reported for each fire pixel. We assume that the pixel 
with the maximum MODIS FRP value ever detected within 
each biome burns completely (i.e., active fire area for that 
pixel equaled 1 km 2 ). We then scale the rest of the MODIS 
FRP values by this maximum and sum those resulting frac- 
tional areas to obtain total active fire area for each plume. 
For example, the maximum FRP value detected over the 
boreal forest is 1710 W/m 2 . Thus, an FRP of 800 W/m 2 
within the same type of vegetation is assigned an active fire 
area of 0.47 km 2 . Finally, the Dual-Band approach esti- 
mates the active fire area using the bi-spectral method of 
Dozier [1981] applied to the individual fire clusters. 

As shown in Table 3, these four approaches provide the 
model with active fire areas that differ significantly, with a 
range factor of 50-100. For example, the WRF-Pxl is about 
10 times the FLAMBE-Pxl. The active fire area mean from 
scaled-FRP is close to the FLAMBE mean, and the Dual- 
Band active fire area is 7-10 times smaller than those. 

3. Results 

3.1. Evaluation of fire input configurations 

Figure 2 presents the relationship between the plume-rise 
model injection heights and the MISR plume heights for the 
various fire input parameters (Table 2 and Table 3). For 
clarity, they are also shown for 500-m bins of MISR plume 
heights. 

It is apparent that the plume model gives different in- 
jection heights depending on the input used, and we ex- 
plore the sensitivity to key parameters below. Overall, 
the plume-rise model does not reproduce well the smoke 
heights observed by MISR as the maximum r 2 obtained is 
0.33. The FLAMBE-Pxl/FRPxlO (Figure 2h) and scaled- 
FRP/FRPxlO (Figure 2p) show the best fit to the MISR 
observations. The WRF-Pxl/FT06 (Figure 2a) and WRF- 
Pxl/FUEL (Figure 2b) configurations give the largest injec- 
tion heights, in particular for the lowest plumes observed 
by MISR. The overestimate for these configurations results 
from the combination of larger fire areas and total heat 
flux values used as inputs. The fire areas for WRF-Pxl are 
more than one order of magnitude larger than those of the 
other three approaches (Table 3). The FT06 and FUEL 
heat fluxes are on average twice as large as the FRPxlO and 
Dual-Band ones, for forest fires, which represent 35% of the 
plumes in the dataset. Based on these results, we discuss 
further the scaled-FRP/FRPxlO (Figure 2p) configuration; 
we consider it the best approach for estimating the fire areas 
and heat fluxes for use in the plume rise model, as it gives 
the best fit to the MISR observations (r 2 =0.33). Addition- 
ally, we discuss further WRF-Pxl/FT06 since it has been 
adopted by the WRF-Chem model framework [ e.g . Pfister 
et al., 2011; Sessions et al., 2011]. 

To better interpret the model results relative to the fire in- 
puts, we examine in more detail the model injection heights 
as a function of active fire area and total heat flux for four of 


the 16 configurations: WRF-Pxl/FT06, WRF-Pxl/FRPxlO, 
scaled-FRP/FT06 and scaled-FRP /FRPxlO, and compare 
them to the observed MISR heights (Figure 3). We also 
show the relationship between total heat flux and fire area 
for WRF-Pxl/FT06 and WRF-Pxl/FRPxlO in Figure 4. 

As shown in Figure 3, MISR plume heights increase with 
fire area, regardless of the method used to estimate the ac- 
tive fire area, although significant variability is present in the 
distributions. The median height of large fires (i.e., > 850 ha 
in WRF-Pxl and > 150 ha in scaled-FRP) are 1.5-1. 8 km 
higher than those of small fires (i.e., < 100 ha in WRF-Pxl 
and < 2.5 ha in scaled-FRP). This is consistent with the 
results of Val Martin et al. [2010], who showed that MISR 
plume heights were related to MODIS FRP with significant 
variability in a log-log relationship, and the fact that we use 
MODIS FRP and fire pixels as the basis for estimating fire 
area. Obviously, similar results are reflected in the model 
injection heights, as the model is sensitive to the fire area 
and heat flux values. For example, WRF-Pxl/FT06 (Fig- 
ure 3a) provides the plume-rise model with large values of 
fire area and heat flux, and the medians of the simulated 
injection heights are typically at least 1 km higher than 
the observations. Conversely, scaled-FRP/FT06 initializes 
the model with lower fire areas than WRF-Pxl/FT06, and 
the simulated median injection heights are on an average 
0.5 km lower than those simulated in WRF-Pxl/FT06, al- 
though still 0.5 km higher than the observations. 

The WRF-Pxl/FRPxlO and scaled-FRP/FRPxlO config- 
urations use lower heat flux values than WRF-Pxl/FT06 and 
scaled-FRP/FT06. In FRPxlO 80% of the cases have heat 
flux values <10 kW/m 2 , compared to only 6% in FT06. 
Figure 4 shows that for WRF-Pxl/FRPxlO a good linear 
correlation (r 2 =0.7) exits between fire area and heat flux 
since both parameters are estimated using FRP values at 
each MODIS fire pixel. Hence, smaller fires (i.e., with fewer 
number of fire pixels and/or lower FRP values) tend to have 
lower fire energy emission rates. Similar relationships are 
found for scaled-FRP/FT06 and scaled-FRP /FRPxlO (not 
shown), with lower values for active fire areas and total 
heat fluxes (Table 3). In the case of WRF-Pxl/FT06 (and 
scaled-FRP/FT06), no relationship exists between fire area 
and heat flux, as only three values, which depend on veg- 
etation type, are used for total heat flux, as described in 
section 2.4.1. 

An important property of the plume-rise model is its abil- 
ity to predict which plumes reach the FT. Figure 5 shows 
the relationship between injection height predicted by the 
model and observed by MISR, and highlights the plumes 
that are located in the FT, for the WRF-Pxl/FT06 and 
scaled-FRP/FRPxlO configurations. Following Kahn et al. 
[2008], we consider a smoke plume to be in the free tro- 
posphere with high confidence when the (plume-BL) height 
difference is > 500 m, to account for errors in both the MISR 
heights and the GEOS BL heights. We also report below 
results when the (plume-BL) height difference exceeds zero. 
The WRF-Pxl/FT06 case implies that 77-94% of the plumes 
reaching the FT (defined as (plume-BL) height>500 m and 
(plume-BL) height>0 m, respectively), whereas the frac- 
tions for the MISR data are 24-48%. As expected, the 
scaled-FRP/FRPxlO configuration implies a lower fraction 
of plumes in the FT, 16-35%, closer to the MISR results. 
However, the model does not put the correct plumes in the 
FT. As shown in Figure 5b, the MISR data implies that 
137 plumes are in the FT, whereas the model puts only 75 
plumes there, and only 42 plumes are in the FT according 
to both MISR and the model. In practice, the model locates 
44% (33 out of 137) of the plumes in the FT and misses 70% 
(95 out 137) of the MISR plumes in the BL 



X - 6 


VALMARTIN ET AL.: SMOKE PLUME-RISE CONSTRAINTS 


3.2. Effect of atmospheric stability on simulated 
injection height 

Val Martin et al. [2010] demonstrated the important ef- 
fect of the atmospheric structure on the ultimate rise of fire 
emissions, as proposed by Kahn et al. [2007]. We illustrate 
the model sensitivity to the heat flux and fire area for two 
very different cases of atmospheric stability for the two rep- 
resentative plumes shown in Figure 1. We consider atmo- 
spheric stability as the change in potential temperature with 
height [Holton, 1992] at the location and time of the MISR 
plumes, as in the work of Kahn et al. [2007] and Val Martin 
et al. [2010]. Figure 6 depicts the atmospheric conditions for 
the two MISR plumes used as examples, and the results of 
the final injection heights as a function of the active fire area 
and total heat flux. We show model results for total heat 
flux ranging from 1 to 100 kW/m 2 and active fire area rang- 
ing from 1 to 250 ha, which are the range of values obtained 
with the scaled-FRP/FRPxlO approach. 

The profile in Figure 6a shows the existence of a weak, 
broad stable layer, with maximum of 5 K/km centered at 
about 3 km altitude and 3 km deep, whereas the profile 
in Figure 6c shows a strong stable layer, with maximum 
of 16 K/Km, centered at 750 m and 1 km deep. These two 
cases exhibit similar wind speed conditions (not shown). For 
the “strong stability” case (Figures 6c, d), the model simu- 
lates a plume height that varies from 0.5 to 2.4 km over the 
range of input studied; for the “weak stability” case (Fig- 
ures 6a, b), the range is from 0.6 to 4.3 km. The model gives 
different injection heights due to the presence of the strong, 
low altitude stable layer in the “strong stability” case that 
causes the simulated plume to lose vertical velocity faster 
than in the “weak stability” case, as explained by Freitas 
et al. [2007] in a sensitivity test performed with two differ- 
ent thermodynamic conditions. In that work, Freitas et al. 
[2007] also highlighted the weak sensitivity of the model to 
the heat flux values. We find, however, that the model also 
has a low sensitivity to fire area because a given plume height 
is associated with multiple combinations of these two input 
parameters. For example, in the “weak stability” case, a 
plume height of 3 km can be obtained with model simula- 
tions that use fire area ranging from 25 to 150 ha, and total 
heat flux from 25 to 100 kW/m 2 . 

Based on a simple 1-D convective parcel model, Kahn 
et al. [2007] showed that, in addition to atmospheric stabil- 
ity structure and total heat flux, entrainment is a complex 
parameter associated with atmospheric stability that can af- 
fect significantly the final rise of the plume. As described in 
Freitas et al. [2010] and mentioned in section 2.2, two en- 
trainment parameters are included in the plume rise model: 
lateral entrainment, to account for the loss of buoyancy due 
to entrainment of cold air into the plume; and dynamic en- 
trainment, to account for the horizontal displacement of the 
plume due to wind drag. Both parameters are inversely pro- 
portional to the plume radius. In the plume-rise model, the 
formulation for the lateral entrainment coefficient is 2oR _ , 
in which R is the radius of the plume and a is assumed equal 
to 0.05 [ Freitas et al, 2010]. We tested slightly smaller and 
larger entrainment coefficients (i.e., a of 0.025 and 0.1) us- 
ing the scaled-FRP/FRPxlO configuration. The mean injec- 
tion heights were 14% larger with a equal to 0.025 and 12% 
smaller with a equal to 0.1. However, the change in entrain- 
ment coefficients did not significantly improve the correla- 
tion between the model injection heights and the observed 
MISR plume heights (that is, the r 2 was about 0.3). 

3.3. Constraints for smoke plume rise 

As shown in section 3.1, constraining quantitatively the 
plume-rise model with observations is a difficult task due to 
the uncertainty in the input parameters. Clearly, some pre- 
vious assumptions adopted in 3-D models do not match the 
MISR plume heights at all well. Here we propose a simpli- 
fied approach for computing plume injection heights using 


two main physical factors: buoyancy flux and atmospheric 
stability structure. This analysis is based on the findings 
presented in Val Martin et al. [2010] for the entire MISR 
dataset. The authors showed that plumes reaching the FT 
have larger FRP values that those trapped within the BL, 
and that the depth of the BL and the presence of stable 
layer in the FT influence the smoke plume height. To illus- 
trate our simple approach, we classified the MISR plumes 
by altitude and analyzed the atmospheric stability struc- 
ture and fire radiative conditions within this classification. 
Figure 7 displays the results for three cases: observations 
of low-altitude (<1000 m), medium-altitude (1000-2500 m) 
and high-altitude (>2500 m) plumes, and Table 4 summa- 
rizes the statistics of the parameters in these categories for 
all observations and by biome type. This analysis shows 
that smoke plumes reaching high altitudes exhibit statisti- 
cally higher FRP than those at low altitude, with a mean of 
about 2100 MW higher. In addition, plumes reaching high 
altitudes are associated with weaker atmospheric stability 
conditions, and plumes with low altitude get confined below 
the BL, consistent with Val Martin et al. [2010]. For ex- 
ample, the atmospheric stability of high plumes (Figure 7f) 
is most often lower than 10 K/Km, whereas it frequently 
reaches more than 10 K/km in low-altitude plumes (Fig- 
ure 7b). In addition, a distinct stable layer at around 1500 m 
is typically present in low plumes, and smoke never punches 
through this layer, which represents the BL cap. In contrast, 
a distinct layer, if present, appears at about 2000-2500 m 
in high plumes, and smoke often punches through. Results 
from medium-altitude plumes (Figure 7d) show intermediate 
behavior. We replicated these analyses stratified by biome, 
and found similar patterns, with different FRP values, rep- 
resentative of forest and shrubland fires (large FRP) and 
grassland fires (small FRP). However, further subdivision of 
the data reduced the sample sizes, and some results were 
not statistically robust, in particular for grassland fires (Ta- 
ble 4). 

In Figure 8, we classify the observations further by at- 
mospheric stability and plume height, to investigate the fire 
radiative heat under these conditions. We consider highly 
stable those atmospheric profiles having stable layers with 
stability above 10 K/km (e.g, Figure 6c), weakly stable 
those atmospheres with stability below 10 K/km (e.g., Fig- 
ure 6a), and neutrally stable those without a distinct sta- 
ble layer. Figure 8 shows that high FRP values are as- 
sociated with stronger atmospheric stability in each plume 
elevation category. Fires with low altitude smoke are as- 
sociated with larger FRP values in highly stable conditions 
(246 MW) than in weakly stable (164 MW) and neutrally 
stable (136 MW) ones. A similar pattern is found for fires 
with smoke that reaches altitudes of 1000-2500 m (867, 412 
and 316 MW, respectively). We have fewer observations for 
fires with high altitude smoke, so the distributions are not 
statistically representative. However, it is clear in Figure 8 
that the more stable the atmosphere, the larger the FRP 
associated with a given height category, suggesting greater 
buoyancy flux is required for the smoke to reach similar al- 
titudes, as expected on physical grounds. 

4. Discussion and Conclusions 

We used MISR plume heights over North America to eval- 
uate the ability of the plume-rise model described in Freitas 
et al., (2007,2010) to produce accurate injection heights. We 
initialized the model with GEOS assimilated meteorological 
fields and estimated fuel moisture content at the locations 
and times of the MISR observations. We tested the model 
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with four estimates each for active fire area and total sensi- 
ble heat flux. The analysis shows that the plume-rise model 
gives different injection heights depending on the input used, 
and is not able to reproduce the plume heights observed by 
MISR uniformly, over the range of conditions studied (max- 
imum r 2 is 0.33 in Figure 2). In addition, the model fails to 
locate the right plumes above the BL (Figure 5), although 
this is crucial information needed for atmospheric models to 
simulate properly the dispersion of fire emissions. 

Val Martin et al. [2010] analyzed a 5-year record of MISR 
plume heights over North America (3367 plumes) using dif- 
ferent definitions of plume heights. Our analysis showed 
that, using the best estimated median of the vertical distri- 
bution of stereo-height retrievals for the entire plume, 4-12% 
of the plumes reach the FT ( i.e ., plume-BL height>500 m), 
depending on the year. These fractions were higher, 15- 
28%, for the best estimated maximum plume height of the 
distribution. Using a less conservative definition for a plume 
residing in the FT (i.e., plume-BL height>0 m), the frac- 
tions were 17-30% and 38-55% for the best estimated me- 
dian and maximum height, respectively. Here, we modified 
our previous definition of plume height to be compatible 
to the injection height simulated by the plume-rise model, 
and considered the best estimated maximum height above 
the fire source. Using this plume height definition and a 
smaller number of observations (584 plumes), we find that 
24-48% of the plumes reach the FT (defined as (plume-PBL) 
height>500 m and (plume-PBL) height>0 m, respectively). 
It is important to note that the maximum plume height defi- 
nition does not imply that the entire plume is above the BL, 
as plumes are typically about 850 m thick and some fraction 
of the smoke may remain in the BL [ Val Martin et al., 2010]. 

Our plume-rise model simulations using the scaled- 
FRP/FRPxlO configuration show that these input settings 
initially provided the best fit to the observed MISR plume 
heights, since they give a fraction of plumes in the FT closest 
to that observed by MISR (16-35% versus 24-48%). How- 
ever, a more detailed analysis (Figure 5b) indicates that 
the model incorrectly locates 44% of the plumes in the FT, 
whereas it misses 70% of the actual cases where the plume 
is in the FT. 

Brioude et al. [2009] used the plume-rise module with 
the FLEXPART Lagrangian Particle Dispersion Model to 
study the impact of biomass burning on clouds off the Cal- 
ifornia coast. In that work, the fire size and heat flux for 
each fire were retrieved from WF-ABBA; for the heat flux, 
they used our approach, which assumes that radiant energy 
is approximately 10% of the total heat energy released by 
the fire. Using these fire settings, and comparing the injec- 
tion heights to the FLEXPART BL, they reported that 30% 
of the plumes reach the FT. However, they did not have 
observations to evaluate the modeled injection heights. 

The current configuration in the WRF-Chem plume-rise 
module (i.e., WRF-Pxl/FT06) does not give the model the 
correct fire inputs. The model overestimates the injection 
height significantly, and locates 77-94% of the plumes in the 
FT (Figure 5a), as determined by the GEOS BL heights. 
This setting uses only five heat flux values taken from con- 
ditions appropriate for Brazil and applies them to fires in 
North America. Sessions et al. [2011] compared injection 
heights simulated by WRF-Chem with an embedded plume- 
rise model to MISR heights obtained during the ARCTAS 
campaign over Siberia and Canada in 2008. They used 
two preprocessing configurations to determine fire locations 
and fire size: one is the default in WRF-Chem and re- 
lies on WF-ABBA and MODIS pixels, the other is based 
on the FLAMBE dataset. In a 10-day run, the WRF- 
Chem placed more than 75% of the plumes above the BL, 
for both configurations. Similarly, Pfister et al. [2011] re- 
ported that the WRF-Chem plume-rise module released 
about 50% of the emissions into the FT in a study over Cal- 
ifornia during ARCTAS-CARB. Using the same definition 


for plumes above the BL (i.e., [plume height-BL]>0 km), our 
analysis indicates that a larger fraction of plumes reaches 
the FT (94%). As discussed in detail by Sessions et al. 
[2011], WRF-Chem BL heights are lower than GEOS BL 
heights due to different BL schemes. In addition, analysis 
of plumes from different locations, different plume model 
versions (with and without wind effects) and different mete- 
orological fields may be also responsible for the differences 
in our results regarding the fractions above the BL. In any 
case, it is clear that the plume-rise module package sup- 
ported by WRF-Chem overestimates the injection height 
compared with MISR observations, that indicate typically 
about 24-48% of plumes reach the FT over North America. 

Based on previous studies [ e.g . Kahn et al., 2007; Freitas 
et al, 2007], it is clear that smoke plume rise is governed by 
buoyant energy produced by the fire, ambient atmospheric 
stability and entrainment. In addition, latent heat release 
and wind are additional factors that can affect plume rise 
[ Freitas et al, 2010]. However, our analyses show that con- 
straining the plume-rise model quantitatively is a difficult 
task, due to uncertainty in the input parameters. Since 
the buoyant energy varies significantly depending on the fire 
type, and typically exhibits strong heterogeneity across ac- 
tive burning areas, it is difficult to measure directly. The 
model assumes that sensible heat flux is 55% of the total 
heat flux, and we further infer total heat flux from instan- 
taneous fire radiative measurements made by MODIS, con- 
sidering that 10% of the total heat flux is radiant energy. 
However, MODIS FRP measurements may be uncertain for 
several reasons. For example, the presence of dense smoke 
can lead to underestimates of the FRP, and the relationship 
between observed spectral signals and FRP can depend on 
the spectral emissivity assumed in the retrieval, which likely 
varies depending upon the fire dynamics (e.g., smoldering 
versus flaming) and flame thickness. [ Kahn et al, 2007]. 

The atmospheric conditions used to constrain the plume- 
rise model are often obtained from global or regional atmo- 
sphere models, and the coarse spatial and temporal resolu- 
tion can introduce uncertainties in the atmospheric stabil- 
ity calculation and the boundary layer heights [Val Martin 
et al., 2010; Sessions et al, 2011]. Entrainment is a yet 
more difficult parameter to constrain given current knowl- 
edge, as it depends on small-scale interactions between the 
plume and the surrounding environment (Kahn et al., 2007]. 
In the plume-rise model, lateral and dynamic entrainment 
are parameterized with simple coefficients inversely propor- 
tional to the cross-sectional basal area of the plume, which is 
taken from the active fire area estimate. However, fires typi- 
cally have non-circular cross-sections, and burning generally 
does not occur at a uniform temperature. 

Other challenges in modeling fire emissions are the di- 
urnal and seasonal fire intensity variations, which are not 
taken into account by the plume-rise module used in this 
study, as it relies uniquely on five heat flux values depend- 
ing only on vegetation type. Evaluation of the plume-rise 
model is limited in our case to the MISR observation time, 
which are 11:00-13:45 local time for our dataset [Val Mar- 
tin et al., 2010]. The MISR local observing time does not 
coincide with the maximum in fire intensity, which typically 
occurs in the afternoon (e.g. Giglio et al., 2006; Vermote 
et al, 2009; Mu et al, 2011]. The MODIS FRP diurnal 
cycle indicates that fire intensity peaks between 13:00 to 
17:00 in a variety of biomes [Giglio, 2007; Roberts et al, 
2009; Vermote et al, 2009]. Furthermore, Val Martin et al. 
[2010] compared the variability of MISR plume heights and 
MODIS FRP on a monthly basis, and concluded that the 
annual variability of plume heights was driven mainly by the 
annual cycle of fire radiative heat. 
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The use of a plume-rise model requires the calculation of 
an injection height for each fire, obtained by numerical inte- 
gration of a system of equations. This numerical approach 
makes the process computationally expensive [Freitas et al., 
2007; Guan et al., 2008]. In addition, the validation of the 
plume-rise model has been based to date on a few case stud- 
ies, with the exception of the work of Sessions et al. [2011]. 
Here, using a large plume dataset and several approaches 
to estimate fire properties, we demonstrate that the model 
currently does not simulate properly the plume heights, and 
this implies that embedding the plume-rise model using the 
currently available fire constraints in large-scale atmospheric 
models is not a realistic approach. 

Our analysis of the main physical constraints for smoke 
plume rise (Figures 7 and 8, and Table 4) shows that smoke 
plumes reaching high altitudes exhibit higher FRP and 
weaker atmospheric stability conditions than those at low 
altitude, which tend to be confined below the BL. The sim- 
plified relationship of FRP and atmospheric stability, sum- 
marized in Table 4, can be used as a basis to develop an inex- 
pensive, empirically based parameterization for computing 
the plume injection heights for fires. It is important to note 
that entrainment is also a key factor that drives plume rise, 
but we could not constrain it with the information available 
to date. 

A number of FRP-based emission inventories have been 
developed recently, e.g., the Global Fire Assimilation Sys- 
tem (GFAS) [Kaiser et al, 2012] and the Quick Fire Emis- 
sion Dataset (QFED) [da Silva, NASA, personal communi- 
cation], and successfully incorporated into CTMs. For ex- 
ample, QFED is operational in the GEOS-5 aerosol fore- 
casting system and Goddard Chemistry Aerosol Radia- 
tion and Transport (GOCART) model [Petrenko et al., 
2012], whereas GFAS is used within the operational fore- 
casting system being developed for European monitoring of 
the global atmosphere (http://www.gmes-atmosphere.eu/). 
Most commonly, however, CTMs use biomass burning emis- 
sion inventories derived from satellite measurements of post- 
fire burned area and estimates of biomass consumed per 
unit area, e.g., the Global Fire Emissions Database (GFED) 
[van der Werf et al, 2010]. A number of studies have shown 
a direct linear relationship between FRP and biomass con- 
sumed [e.g. Wooster et al., 2005; Freeborn et al, 2008], and 
a set of empirically derived factors have been developed to 
relate satellite-derived FRP to combustion rates from dif- 
ferent biomes [e.g. Kaiser et al, 2012]. These factors may 
be used inversely to connect the fuel consumption values re- 
ported in the GFED inventory to FRP values, on a monthly 
or daily basis. Thus, CTMs could use the information on 
atmospheric stability and FRP or a direct relationship to 
FRP to treat empirically the injection heights on large-scale 
studies. 

Another simplified parameterization yet to be tested is 
based on the statistical analysis of MISR plume climatolo- 
gies for extratropical [Val Martin et al, 2010] and tropical 
[Tosca et al., 2011] biomes. Based on the MISR observa- 
tions, we recommend injecting about 15% of the emissions 
into the FT for extra-tropical forest and shrubland fires and 
about 5% for tropical forest and grassland fires. Fire emis- 
sions released in the BL will be mixed throughout the BL, 
whereas emissions released in the FT will be injected at the 
height of the stable layer or spread-out vertically in those 
model vertical grid boxes without distinct atmospheric sta- 
ble layers, following the work of Val Martin et al. [2010]. 

To reduce the uncertainty associated with plume rise 
models, we recommend small-scale, simultaneous field mea- 
surements of all the parameters driving smoke plume-rise. 
The necessity of new field campaigns to quantitatively char- 
acterize plume rise has also been recommended by Ichoku 
et al. [2012] after a comprehensive review of satellite obser- 
vations on biomass burning for atmospheric modeling. In 
future work, we will explore the implementation of the sim- 
plified injection height parameterizations within a CTM. 


Appendix A: Development of a fuel consumption 
map over the North America boreal region 

Previous studies have shown that current estimates for 
fuel consumption in boreal forest fires are underestimated 
[e.g. Amiro et al., 2001; Kasischke and Bruhwiller, 2002; 
Turquety et al., 2007]. We developed a new map of fuel con- 
sumption using the U.S. Forest Service (USFS) CONSUME 
model (version 3) [Ottmar, 2009], adapted for boreal ecosys- 
tems. This model uses as input fuel loading and moisture 
conditions, and produces estimates of fuel consumption. It is 
designed to use the fuel loadings for the Fuel Characteristic 
Classification System (FCCS) for the U.S. (lower 48 states 
and Alaska) [Ottmar and Baker, 2007]. CONSUME requires 
maps of fuelbeds. For Canada, fuelbed information at 1-km 
resolution was obtained from the Canadian Fire Behavior 
Prediction (FBP) system, which is based on satellite-derived 
vegetation data [Nadeau et al, 2005]. For Alaska, we used 
a map created by the USFS, which follows the classification 
scheme of Nadeau et al. [2005] [Martin Parker, USFS, per- 
sonal communication]. We mapped the Canadian fuelbeds 
onto their corresponding type in the U.S. FCCS [Ottmar 
and Baker, 2007] . Some assumptions were needed to match 
some of the Canadian FBP fuelbed types into the U.S. FCCS 
classification system [Nancy French, MTRI, personal com- 
munication], and we summarize those in Table 5. 

To compute fuel consumption, we then imported the 
FCCS fuel loading values for each fuelbed type into CON- 
SUME. The effect of burning duff layer was considered 
within the fuel consumption values in the fuelbed, as CON- 
SUME takes into account burning of the organic soil layer. 

We adjusted the organic soil layer in some of the fuelbed 
types for that purpose. For example, we increased the 
duff layer depth for Black spruce-feathermoss (FCCS 87) in 
CONSUME from 15 cm to 30 cm to match the description 
of deep, compacted organic layer from the Canadian FBP 
system (FBP C2). 

To represent the changes in fuel moisture regime through- 
out the fire burning season, we computed fuel consump- 
tion values in each fuelbed type for five different types 
of moisture conditions ( i.e ., wet, moist, moderately-dry, 
dry and extra-dry) and considered changes in moisture 
conditions in the three fuelbed stratum given in CON- 
SUME (he., woody fuels of 0.6-3 cm diameter (10-hr), 
woody fuels of 8-20 cm diameter (1000-hr) and duff). Ta- 
ble 6 summarizes the fuel moisture percentages for the 
fuelbed stratum in each moisture category. These fuel 
moisture conditions are based on typical conditions used 
in CONSUME and the Fire Emission Production Simula- 
tor (http://www.fs.fed.us/pnw/fera/feps/index.shtml); val- 
ues are consistent with previous reports of soil moisture con- 
tent, based on field measurements over the boreal region 
[Ottmar and Baker, 2007; Jandt et al, 2005]. 

We derived fuel consumption maps at 1° latitude by 
1° longitude for the five fuel moisture conditions, for use in 
climate and atmospheric models. Figure 9 shows the boreal 
region fuel consumption map for moderately-dry conditions 
as an example. We also include in the map the fuel consump- 
tion values for the lower 48 states described in Spracklen 
et al. [2009]. Table 3 summarizes the fuel consumption val- 
ues for the fuelbed types by moisture conditions. In this 
work, to determine the total heat flux values (in W/m 2 ), we 
use the fuel consumption values obtained for moderately-dry 
conditions (Figure 9) and a standard heat of combustion for 
extra-tropical vegetation of 18.7e6 J/Kg [Trentmann et al., 

2006], and assuming that the average burning time of the 
fire is two hours. 
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Figure 1. Example of the vertical distribution of MISR stereo-height pixels for a boreal 
forest fire in central Canada (a) and a grassland fire in Texas (c). Also shown are 1-D 
plume-rise model outputs obtained for the forest fire with heat flux of 120 kW/nt 2 and 
area of 127 ha (b), and for the grassland fire with heat flux of 55 kW/m 2 and fire area of 
3.5 ha (c). 
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Figure 2. Relationship between MISR observed maximum plume heights and 1-D 
plume-rise model simulated plume heights for the active fire area configurations (a)-(d) 
WRF-Pxl, (e)-(h) FLAMBE-Pxl, (i)-(l) Dual-Band and (m)-(p) scaled-FRP, and the total 
heat flux configurations (a),(e), (i) and (m) FT06, (b), (f), (j), (n) FUEL, (c), (g), (k), 
(o) Dual-Band and (d), (h), (1), (p) FRPxlO. See table 2 and text for explanation. 
Regression lines for correlations with r 2 below 0.25 are not shown. Grey lines are the 
two-sided regression line [Ayers, 2001]; black line is the 1:1 relationship. Bar plots are 
the distributions of model heights and 500-m resolution MISR heights. The medians (red 
circles) are shown along with the central 67% (box) and the central 90% (thin black lines). 
Red lines are the regression line of the distribution medians. 
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Figure 3. Distribution of the MISR plume height (green) and the 1-D plume-rise 
model height for (a) WRF-Pxl/FT06 (blue) and WRF-Pxl/FRPxlO (light blue) and (b) 
scaled-FRP/FT06 (blue) and scaled-FRP/FRPxlO (light blue) configurations grouped by 
active fire area. The medians (red circles) and the means (black squares) are shown along 
with the central 67% (blue, light blue or green) and the central 90% (thin black lines). 
The number of observations (in black) and the median heights (in colors) included in each 
distribution are given at the top of the plot. 
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Active Fire Area (ha) 

Figure 4. Relationship between active fire area and total heat flux for WRF-Pxl/FT06 
(in red) and WRF-Pxl/FRPxlO (in black). Black line is the two-sided regression line for 
WRF-Pxl/FRPxlO. 
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Figure 5. Relationship between MISR plume height and the plume-rise model height 


for WRF-Pxl/FT06 (a) and scaled-FRP/FRPxlO (b) configurations. All observations are 


shown in grey, plumes located in the FT by the model are shown in red, and plumes in 


the FT by MISR are shown with open black circles; black line is the 1:1 relationship. 
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Figure 6. Atmospheric stability profiles for the MISR plumes shown in Figure la (a) 
and Figure lb (c), and the associated contour plots for the injection height as a function 
of fire area and total heat flux simulated by the model (b and d) 
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Figure 7. Vertical distribution of the MISR stereo-height retrievals for the 584 plumes 
analyzed, classified as low altitude plumes (a), medium altitude plumes (c) and high 
altitude plumes (e). The atmospheric stability structures associated with this classification 
(b, d, f) are shown as the median (solid line) and the 5th and 90th percentiles (dotted 
lines). Numerals in red and black indicate the mean of the MODIS FRP and the number 
of plumes in each distribution, respectively. 
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Figure 8. Distribution of MODIS FRP by atmospheric stability (Hg=highly, 

Wk=weakly and Nt =neutrally) for plumes below 1000 m (yellow), plumes between 1000- 
2500 m (blue) and plumes above 2500 m (green). (See text for further explanation). The 
medians (red circles) and the means (black squares) are shown along with the central 67% 
(in color) and the central 90% (thin black lines). The number of observations (in black) 
and the median heights (in colors) included in each distribution are given at the top of 
the plot. A distribution with 2 plumes is shown as the average (black square). 
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Figure 9. Fuel consumption (in kg dry matter/m 2 ) map for North America for 
moderately-dry fuel moisture conditions in the boreal region. 
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Table 1. 1-D Plume- rise models used for smoke injection height simulations. 


Plume-rise Parameterization 

Scheme 

Hosted Model 

Briggs [1969] 

Empirical 

CMAQ a 

HYSPLIT 6 

Lavoue et al. [2000] 
Kahn et al. [2007] 

Achtemeier et al. [2011] 
Freitas et al. [2007, 2010] 

Empirical 

Physical 

(diagnostic) 

Empirical 

Physical 

(prognostic) 

CMAQ C 

CATT-BRAMS d 

NCAR-CAM 6 

FLEXPART^ 



WRF-Chem 9 

Rio et al. [2010] 
Sofiev et al. [2012] 

Physical 

(prognostic) 

Empirical 

LMDZ 


a Raffuse et al. [2012]; b Stein et al. [2009]; c Achtemeier et al. [2011] and 

Liu et al. [2010] \ d Freitas et al. [2007]; e Guan et al. [2008 }■/ Brioude et al. [2009]; 

9 Pfister et al. [2011], Grell et al. [2011] and Sessions et al. [2011] 
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Table 3. Summary of the active fire area and total heat flux for each biome a . 


Configuration 

Forest 

Biome Type 
Shrubland 

Grassland-Crop 



Fire Size (ha) 


WRF-Pxl 

445 ± 467 (250) 

410 ± 375 (250) 

230 ± 208 (150) 

FLAMBE-Pxl 

45 ± 47 (25) 

41 ± 38 (25) 

23 ± 21 (15) 

scaled-FRP 

61 ± 101 (20) 

44 ± 81 (18) 

20 ± 25 (8) 

Dual-Band 

8 ± 15 (3) 

6 ± 11 (3) 

2 ± 3 (15) 


Heat Flux (kW/m 2 ) 


FT06 6 

30-80 

4.4-23 

3.3 

FUEL 

48 ± 30 (52) 

36 ± 24 (38) 

32 ± 27 (40) 

FRPxlO 

10 ± 16 (3) 

7 ± 13 (3) 

3 ± 4 (2) 

Dual-Band 

23 ± 19 (18) 

18 ± 14 (14) 

20 ± 19 (14) 


“Reported are mean, standard deviation (SD) and median, as MeaniSD (Median). 
6 Shown upper and lower limit values from Freitas et al. [2006]. 
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Table 4. Statistical summary of the parameters shown in Figure 7 for all observations and 


for each biome“. 


Plume 

Height 6 

FRP (MW) 
Mean±SD Median 

Stability 

Stable layer (K/km) 

Altitude (m) 

N 

Low 

318 ±485 

164 

All Observations 
5.3-13.2 

1800 

225 

Medium 

785 ±997 

409 

5.4-14.0 

1800 

286 

High 

2478±3070 

1450 

5. 1-9.4 

2800 

73 

Low 

428 ±672 

226 

Forest 

4.8-13.6 

1200 

81 

Medium 

948 ±1230 

409 

5.5-14.7 

1800 

86 

High 

3032±3130 

2170 

4.8-12.1 

2300 

29 

Low 

259 ±326 

147 

Shurbland 

5.7-12.6 

1800 

126 

Medium 

744 ±899 

425 

5.3-11.4 

1800 

182 

High 

2250±3170 

1370 

5. 2-9.4 

2800 

39 

Low 

177 ±192 

136 

Grassland- Crops 
5.7-13.2 

1700 

17 

Medium 

377 ±437 

239 

5.8-14.6 

1800 

17 

High 

882 ±775 

957 

9.2-13.1 

2800 

3 


“For FRP, mean, standard deviation (SD) and median are reported; for stability, the median 


and the 90th percentile of the absolute stability, and the altitude of the stable layer are reported; 


N is number of observations. Results for each biome with more than 50 data points are shown 


bold-faced. 

6 Plume heights are ’’low” for plnmes<1000 m, ’’medium” for plumes between 1000 and 2500 m 


and ’’high” for plnmes>2500 m (see text for explanation). 
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Table 6. Fuelbed stratum moisture (in %) for each moisture category. 
Moisture Fuelbed Stratum 


Conditions 

10-hr 

1000-hr 

duff 

Extra-Dry 

6 

8 

25 

Dry 

8 

12 

40 

Moderately-Dry 

9 

15 

70 

Moist 

12 

22 

150 

Wet 

20 

31 

250 


Table 7. Summary of the fuel consumption values (in kg/m 2 ) by fuelbed and fuel moisture 


content conditions. 


Fuelbed 

Extra-Dry 

Dry 

Moderately-Dry 

Moist 

Wet 

Cl 

6.73 

6.01 

4.40 

1.19 

1.08 

C2 

15.15 

14.02 

11.53 

4.41 

2.18 

C3 

15.00 

14.95 

14.80 

13.90 

9.31 

C4 

5.61 

5.41 

5.07 

4.70 

4.26 

C5 

14.97 

14.79 

14.41 

12.89 

7.04 

C6 

8.90 

8.27 

7.00 

5.74 

3.56 

C7 

5.21 

4.93 

4.46 

3.87 

3.68 

D1 

8.78 

8.40 

7.30 

4.36 

4.08 

M1-M2 

4.31 

3.80 

3.18 

2.90 

2.77 

01 

2.36 

2.07 

1.44 

1.25 

1.13 

Ola/Olb 

1.25 

0.85 

0.40 

0.40 

0.34 

Tundra 

3.62 

1.65 

1.62 

1.26 

0.91 
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